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GNSS for Precision 
Frequency and Timing 


¢ “Board level” GPS receivers have offered a pulse-per-second 
(“PPS”) output for many years 

¢ Even old receivers can provide accuracy within 100 
nanoseconds; toward 10 ns with a lot of effort 

¢ Useful as reference for GPS disciplined oscillator 
- But GPS PPS is “noisy” in the short term; it has second-to-second jitter 

worse than a local crystal oscillator 

— More on that in my other presentation 


What’s Changed? 


¢ New devices incrementally better 


The big news: 
affordable dual-freq receivers 
- Used to cost $$$$$ 


- Compensate for ionospheric factors 


- Use higher-resolution codes 
—- Allow raw phase output 
ZED-F9 units <$200 

How does this affect timing 
performance? 


L1 CARRIER 1575.42 MHz 
L1 SIGNAL 


C/A CODE 1.023MHz 


NAV/SYSTEM DATA 50 Hz 
C) Modulo 2 Sum 


+4 L2 SIGNAL 
GPS SATELLITE SIGNALS 


P-CODE 10.23 MHz 
— 


L2 CARRIER 1227.6 MHz 


Time to Test! 


e Thanks to the HamSci consortium, | was able to get my hands 


on several 


Model 


LEA-M8F 
NEO-M8N 
NEO-M8P 
NEO-M8T 
NEO-M9N 
ZED-F9P 
ZED-F9OT 


Current u-blox receivers:* 
Features 


L1, frequency and time sync, disciplined oscillator, RAWX, no gErr 
L1, navigation, no qErr, no RAWX, no OD 

L1, positioning, qErr, RAWX, RTK engine 

L1, timing, qgErr, RAWX, no high-precision IIh output 

L1, navigation, no qErr, no RAWX, no OD 

L1/L2, positioning, qErr, RAWX, RTK engine 


L1/L2, timing, qErr, RAWX, no high-precision Ilh output 


*Support to this project from NSF Grants AGS-2002278, AGS-1932997, and AGS-1932972 is gratefully acknowledged. 


The Test Plan 


Compare the PPS outputs (which u-blox calls “TIMEPULSE” 
outputs) from all receivers against atomic clock 

Measure all receivers simultaneously, using the same 
antenna, to allow direct comparison 

e Measure for several days to get long-term data 

¢ Do additional runs targeting specific capabilities/options 


Step 0: Build the Test Equipment 


¢ About six months spent 

designing/building/coding the 

“multi-TICC” to enable the testing 

- 4 TICC* counters in a box, linked to 
a common timebase and computer 
logging system 

- Allows measuring 8 inputs 
simultaneously with 60 ps 
resolution 


- See: 
https://tapr.wpengine.com/tapr-file-area/time- 
freq/multi-TICC App Note _2020-01.pdf 


* httos://tapr.org/product/tapr-ticc/ 


Step 1: Capture Data 


Receivers set to default values except: 

—- Set to “OD” or “Timing” mode where applicable 

- Fixed location set to surveyed antenna position where applicable 
Receivers connected to common antenna through splitters 
TIMEPULSE outputs connected to multi- TICC inputs 

10 MHz from HP 5071A cesium standard to multi-TICC 
reference input 

multi- TICC output logged to computer 

Collected PPS data for just under 6 days (510K samples per 
receiver) and analyzed with “TimeLab” software 


Overview of Results 


¢° Overview shows two distinct groups, but differences are 
hidden in the noise, so let's look at each group separately 


Allan Deviation oy(t) 


Receiver ADEV @ 1 Sec 


LEA-M8F 4.07x10° 
NEO-M8N __9.03x10° 
NEO-M8P__7.71x10° 
NEO-M8T  9.99x10° 
NEO-M9N _ 9.02x10° 
ZED-F9P  3.85x10° 


ZED-F9T 4.13x10° 


ADEV @ 10K Sec 


6.20x10°° 
1.45x10° 
1.12x10°? 
AZ xA0r= 
1.10x10°" 
SLO 


4.97x10°° 


The ZED-F9 Receivers 
(or, “Why We’re Here Today”) 
“P* model has RTK engine; “T” model has extra I/O and a bit 


lower cost (in quantity) 
Half an order of magnitude better than “8” series! 


Allan Deviation oy(T) Original Phase Difference (Linear residual) 
Averaging window: 1 second 
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Questions Answered: 
Does Sawtooth Correction Work? 


Allan Deviation Oy(T) Original Phase Difference (Linear residual) 
Averaging window: 1 second 
1E-7 
+1.00E-8s 
1E-8 +8.00E-9s 
+6.00E-9s 
1E-9 


+4,00E-9s 
1E-10 

+2.00E-9s 

tea = 0.00E+0s 

-2.00E-9s 


1E-12 
-4.00E-9s 

1E-13 
-6.00E-9s 

1E-14 
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[Trace Cd] CNotes | Sample Interval | Acquired | Instrument] ~1,00E-8s 
ZED-FST Raw PPS HP 5071A 740000 pts 
Se Med : 334795 3495s 335128 335295 335465 335625 335795 335965 336138 336295 3646s 
ZED-F9T Corrected PPS vs HP 5071A 140000 pts 
NEO-MBT Corrected PPS vs HP 5071A 140000 pts [Trace | Notes | Sample Interval | Acquired [Instrument _| 


eg ead ZED-FST Raw PPS vs HP 5071A 140000 pts | TCC 
SEW EES we Ens ZED-FST Corrected PPS vs HP 5071A ha 140000 pts__|_TICC 


Questions Answered: 
Is “OD” Important for Timekeeping? 


e Yes 


Allan Deviation oy(T) Original Phase Difference (Linear residual) 
Averaging window: 1 second 
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[trace | Notes | Sample Interval | Acquired | Instrument _| 


[trace | Notes | Sample Interval | Acquired | Instrument _| 


Questions Answered: 


10 MHZ TIMEPULSE as an RF Source? 
e No 


Phase Noise £(f) in dBc/Hz 


mWAY, 
ST LA 
wie ai 


NEO-M8N (1 Hz) 
NEO-M8T (1 HZ) 


NEO-MON (1 Hz) 
NEO-MON (25 Hz) 
ZED-FST (1 Hz) 
ZED-FST (8 Hz) 


10 Hz 


vs. HP 5065A 
vs. HP 5065A 
vs. HP 5065A 
vs. HP 5065A 
vs. HP 5065A 
vs. HP 5065A 


100 Hz 


10.000 MHz 
10.000 MHz 
10.000 MHz 
10.000 MHz 
10.000 MHz 
10.000 MHz 


7200000 pts 
7200000 pts 
7200000 pts 
7200000 pts 
7200000 pts 
7200000 pts 


10 kHz 


TimePod 5330A 
TimePod 5330A 
TimePod 5330A 
TimePod 5330A 
TimePod 5330A 
TimePod 5330A 


Ap 
16 dB’ 


REF 16,8 dBm ATTEN 28 dB 


26 dB 


MK R 


16,6626 MHz 
-.286 dBm 


Tita 


A 


LCOURTIN 


CENTER 16,8860 MHz 
RES BW 


18 kHz 


SPAN 5@@,@ kHz 
SWP 38 msec 


‘ Autonomous Positioning Performance 


12-hour test logging position data from common antenna 
“CEP” = Circular Error, “EP” = Elevation Error Probabilities 


Receiver CEP — 50% CEP —- 95% CEP — 99% EP —- 50% 
LEA-M8F 1.034 2.900 3.735 1.743 
NEO-M8N 1.090 2.302 3.004 leone 
NEO-M8P dar iire Zolg 2.607 1748) 
NEO-M8T 1.264 2.271 3.441 Wale 
NEO-M9N 0.820 1.684 1.944 12S) 
ZED-F9P 0.559 1.449 1.633 O:817 
ZED-F9OT ord) 2.457 3.142 122° 


This ain't right! 


Positioning Performance 


e “P” receivers have built-in RTK correction processing engine 
—- Provide source of correction data, get mm-level data out 

¢ ~8 hour data collection using Ohio DOT reference station 
network for corrections 

° RTK certainly works 
- Why single-freq M8P performance almost equals dual-freq F9P is 

unknown; M8P is performing better than it should! 

—- Consider these results preliminary 


Receiver CEP — 50% CEP —- 95% CEP — 99% EP — 50% 


NEO-M8P 0.010 0.029 0.035 0.021 
ZED-F9P 0.013 0.025 0.033 0.016 


Post-Processing Performance 


° Collected data from NEO-M8P (GPS), ZED-F9P (GPS), 
ZED-F9T (GPS+GLONASS) as well as survey-grade unit 
¢ Sent off for post-processing* 
° Conclusions: 
— M8P gets to around % meter 
- Dual-freq gets to a handful of mm and competes with survey rx 
- Using GPS+GLONASS improves ZED-F9 results vs. GPS only 


NetRS NEO-M&8P ZED-F9P ZED-F9T 


(GPS) (GPS) (GPS) (GPS+GLONASS) 
24 Hour Sigma (95%) 24 Hour Sigma (95%) 24 Hour Sigma (95%) 24 Hour Sigma (95%) 


LAT ITRF2014 39 xx 42.67100 0.0068 39 xx 42.66852 0.3601 | 39 xx 42.67067 0.0090, 39 xx 42.67086 
LON ITRF2014 -84 xx 41.53109 0.0124  -84 xx 41.53533 0.4131 | -84 xx 41.53164 0.0160 | -84 xx 41.53226 
0.0247 0.6522 0.0370 


EL HGT ITRF2014 247.101 


247.1254 247.1548 


* httos://webapp.geod.nrcan.ac.ca/geod/tools-outils/o00.o0hp 
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